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The J-matrix inverse scattering approach can be used as an alternative to a conventional i?-matrix 
in analyzing scattering phase shifts and extracting resonance energies and widths from experimental 
data. A great advantage of the J-matrix is that it provides eigenstates directly related to the ones 
obtained in the shell model in a given model space and with a given value of the oscillator spacing 
HQ. This relationship is of a particular interest in the cases when a many-body system does not 
have a resonant state or the resonance is broad and its energy can differ significantly from the shell 
model eigenstate. We discuss the J-matrix inverse scattering technique, extend it for the case of 
charged colliding particles and apply it to the analysis of not and pa scattering. The results are 
compared with the No-core Shell Model calculations of 5 He and 5 Li. 



I. INTRODUCTION 

The i?-matrix [1] is conventionally used in the analy- 
sis of scattering data, the parameterization of scattering 
phase shifts and the extraction of resonant energies and 
widths from them. The scattering phase shifts can also 
be analyzed in the J-matrix formalism of scattering the- 
ory [2], 

The inverse scattering oscillator-basis J-matrix ap- 
proach was suggested in Ref. [3]. It was further de- 
veloped in Ref. [4] where some useful analytical formulas 
exploited in this paper, were derived. The J-matrix pa- 
rameterization of scattering phase shifts was shown in 
Ref. [4] to be very accurate in describing NN scatter- 
ing data. This parameterization was used to construct 
high-quality non-local J-matrix inverse scattering NN 
potentials JISP6 [5] and JISP16 [6]. 

In what follows, we demonstrate that the J-matrix can 
be used for a high-quality parameterization of scattering 
phase shifts in elastic scattering of nuclear systems using 
na as an example. The resonance parameters, its energy 
and width, can be easily extracted from the J-matrix 
parameterization. 

Resonance energies are conventionally associated with 
eigenstates above reaction thresholds obtained in various 
nuclear structure models, e. g., in the shell model. This is 
well-justified for narrow resonances, however these eigen- 
states can differ significantly from the resonance energies 
in the case of wide enough resonances. The J-matrix pa- 
rameterization naturally provides eigenstates that should 
be obtained in the shell model or any other many-body 
nuclear structure theory based on the oscillator basis ex- 
pansion (e. g., in the resonating group model) to support 
the experimental nucleon-nucleus scattering phase shifts 
in any given model space and with any given oscillator 
spacing Ml. The shell model eigenstates are provided 
by the J-matrix phase shift parameterization not only in 
the case of resonances, narrow and wide ones, but also in 
the case of non-resonant scattering as well, for example, 
in the case of na scattering in the i + partial wave. We 



will explore these correspondences between the J-matrix 
properties and results from nuclear structure calculations 
in some detail below. 

Next, we extend the oscillator-basis J-matrix inverse 
scattering approach of Ref. [4] to the case of charged 
particles using the formalism developed in Ref. [7] . This 
extended formalism is shown to work well in the descrip- 
tion of pa scattering and the extraction of pa resonance 
energies and widths. The shell model eigenstates desired 
for the description of the experimental phase shifts, are 
also provided by the Coulomb-extended J-matrix inverse 
scattering formalism. 

We also carry out No-core Shell Model [8] calculations 
of 5 He and 5 Li nuclei and compare the obtained eigen- 
states with the ones derived from the J-matrix parame- 
terizations of na and pa scattering. 



II. J-MATRIX DIRECT AND INVERSE 
SCATTERING FORMALISM 

The J-matrix formalism [2] utilizes either the oscillator 
basis or the so-called Laguerre basis of a Sturmian type. 
The oscillator basis is of a particular interest for nuclear 
applications. Here we present a sketch of the oscillator- 
basis J-matrix formalism (more details can be found in 
Refs. [2, 7, 9]) and some details of the inverse scattering 
J-matrix approach of Ref. [4] . The extension of J-matrix 
inverse scattering formalism to the case of charged parti- 
cles is suggested in subsection II B while subsection II C 
describes how to relate the J-matrix inverse scattering 
results to those of the shell model. 

A. Scattering of uncharged particles 

Scattering in the partial wave with orbital angular mo- 
mentum I is governed by a radial Schrodinger equation 

H l u l (E ) r)=Eu l (E,r). (1) 
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Here r = |r|, r = ri — r% is the relative coordinate of 
colliding particles and E is the energy of their relative 
motion. Within the J-matrix formalism, the radial wave 
function ui(E,r) is expanded in the oscillator function 
series 



ui 



{E,r) = J2a nl (E)R nl (r), 



(2) 



n=0 



where the oscillator functions 



Rni(r) = (-1)" 



2nl 



r^T{n + I + 3/2) V^o 



TT-T ) ' (X ) • (3) 
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L" (x) is the associated Laguerre polynomial, the oscilla- 
tor radius ro = y^h/mfl, and m = m\m-2l(m\ + 1112) is 
the reduced mass of the particles with masses m\ and m 2 - 
The wave function in the oscillator representation a n i (E) 
is a solution of an infinite set of algebraic equations 



J2 (H l nn ,-6 nn ,E)a nn (E)=0, 



(4) 



n'=0 

where the Hamiltonian matrix elements H l nn i = 
T l nn , + V^ n , , the nonzero kinetic energy matrix elements 

Tl = ^(2n + Z + 3/2), 

T l n+ i.n = Ti n+1 = - ™y/(n+l)(n + l + 3/2), 



n,n+l 



and the potential energy V 1 within the J-matrix formal- 
ism is a finite-rank matrix with elements 



V 1 , 

y nn' 



v l , 

nn' 





if n and n! < Af; . . 
if n or n! > Af. ^ ' 



The potential energy matrix truncation (6) is the only 
approximation of the J-matrix approach. The kinetic 
energy matrix is not truncated, the wave functions are 
eigenvectors of the infinite Hamiltonian matrix H l nn , 
which is a superposition of the truncated potential en- 
ergy matrix V*, and the infinite tridiagonal kinetic en- 
ergy matrix T l nn , . Note that the Hamiltonian matrix, 
i. e. both the kinetic and potential energy matrices, are 
truncated in conventional oscillator-basis approaches like 
the shell model. Hence the J-matrix formalism can be 
used for a natural extension of the shell model. Note also 
that within the inverse scattering J-matrix approach, 
when the potential energy is represented by the finite ma- 
trix (6), one obtains the exact scattering solutions, phase 
shifts and other observables in the continuum spectrum 
(see [4] for more details). 

The phase shift 61 and the ^-matrix are expressed in 
the J-matrix formalism as 

ta s = Sm{E) - W(£) T^ +1 S^ +1j (E) 
dn 1 C m (E)-g MH (E)T l MM+1 CM + iAEY u 



s = ch\E) w£)n^ + icy +M (g) (8) 

~ C$(25) - G^(E)T^ +1 C^l u (E)' 

where Af+1 is the rank of the potential energy matrix (6), 
the kinetic energy matrix elements T l nn , are given by Eqs. 
(5), regular S n i(E) and irregular C n i{E) eigenvectors of 
the infinite kinetic energy matrix are 



.(-) 



Snl(E) 



irron\ 



T(n + / + 3/2) 



q l+1 exp 



(9) 



C nl (E) = 



tt ro n 



T(n + l + 3/2) r(-i + l/2) 



x exp 



)$(-„-/ -1/2, -l + l/2;q 2 ), (10) 



C lu'( E ) = Cm(E) ± iS n i{E), $(a, b; z) is a confluent 
hypergeometric function, the dimcnsionlcss momentum 

1 = V M- The ma trix elements, 



Af 



Gnn>(E) = -J2 



A=0 



MA)(A|n') 
E\ — E 



(11) 



are expressed through the eigenvalues E\ and eigenvec- 
tors (n|A) of the truncated Hamiltonian matrix, i. e. E\ 
and (n|A) are obtained by solving the algebraic problem 



Af 



Y,H l nn '(n'\\)=E x {n\\), n<Af. 



n'=0 



Only one diagonal matrix element Q\tm{E), 



A=0 



W) 2 ) 

E\ — E' 



(12) 



(13) 



is responsible for the phase shifts and the S-matrix. 
The J-matrix wave function is given by Eq. (2) where 

a nl (E) = cos6 l S nl (E) + sin ,5, C nl (E) (14) 

in the 'asymptotic region' of the oscillator model space, 
n > Af. Asymptotic behavior [2, 7, 9] of functions 
y(E,r) and ^{E.r) defined as infinite series, 



y(E,r) = J2S nl (E)R nl (r) = k 3l (kr) 



n=0 



— > - sm(kr-irl/2) (15) 



and 

00 

^(£,r) = ^C ni (£)i? ni (r) 



-kni{kr) 



n=0 



— > - cos(fcr- nil 2) (16) 

r— >oo r 



3 



[here ji{x) and ni(x) are spherical Bessel and Neumann 
functions, and momentum k = q/r ], assures the correct 
asymptotics of the wave function (2) at positive energies 

E, 

m (E, r) — ► k [cos Si ji (kr) — sin Si ni (kr)] 

r— »oo 

— > - sm[kr + 5i-nl/2]. (17) 

r — >oo T 

In the 'interaction region', n < Af, a n i{E) are expressed 
through matrix elements GnAf(E) (see [2, 7, 9] for more 
details). However a limited number of rapidly decreasing 
with r terms with n < Af in expansion (2) does not affect 
asymptotics of the continuum spectrum wave function. 

A similarity between the J-matrix and i?-matrix ap- 
proaches was discussed in detail in Ref. [7]. Note that 
the oscillator function R n i (r) tends to a 5- function in the 
limit of large n [9, 10], 

R n i(r) — V2r r- 3 / 2 S(r - r c J), (18) 



where 



r c J = 2r 0V /n + ;/2 + 3/4 (19) 



is the classical turning point of the harmonic oscillator 
eigenstate described by the function R n i(r). Therefore 
expansion (2) describes the wave function ui(E,r) at 
large distances from the origin in a very simple man- 
ner: each term with large enough n gives the amplitude 
of ui(E, r) at the respective point r = r^. Within the J- 
matrix approach, the oscillator representation wave func- 
tions a n i{E) in the 'asymptotic region' of n > Af and in 
the 'interaction region' of n < Af are matched at n = Af 
[2, 7, 9]. This is equivalent to the _R-matrix matching 
condition at the channel radius r = b — the J-matrix 
formalism reduces to those of the i?-matrix with channel 
radius b = r^j- if Af is asymptotically large. In particular, 
the function Gmm{E) [see (13)] was shown in Ref. [7] 
to be proportional to the P-matrix (that is the inverse 
R- matrix) in the limit of Af — > oo. 

At small enough values of n, oscillator functions R n i(r) 
differ essentially from the ^-function. Therefore the J- 
matrix with realistic values of truncation boundary Af 
differs essentially from the i?-matrix approach with real- 
istic channel radius values b. It appears that the J-matrix 
formalism with its matching condition in the oscillator 
model space, is somewhat better suited to traditional nu- 
clear structure models like the shell model. 

In the inverse scattering J-matrix approach, the phase 
shifts Si are supposed to be known at any energy E and 
we are parameterizing them by Eqs. (7), (9), (10), and 
(13), i. c. one should find the eigenvalues E\ and the 
eigenvector components {Af\\) providing a good descrip- 
tion of the phase shifts. If the set of E\ and (Af\X) val- 
ues is known, i. e. the function GtfAf(E) is completely 
defined, the S'-matrix poles are obtained by solving nu- 
merically an obvious equation, 



where solutions for q (or E = ^q 2 Ml) should be searched 
for in the desired domain of the complex plane. 

Knowing the phase shifts Si in a large enough energy 
interval < E < E max , one gets the set of eigenenergies 
E\, A = 0, 1, ... , Af by solving numerically the equation 



(21) 



where a_sr+i,i(E) is given by Eq. (14). The equation (21) 
has exactly Af + 1 solutions. The last components (AT| A) 
of the eigenvectors (n\X) responsible for the phase shifts 
and the S'-matrix, are obtained as 



\(AT\\)\ 2 = 



a i T Af,Af+l 



where 



dajsf + ij{E) 



dE 



(22) 



(23) 



E=E X 



The physical meaning of the Eqs. (21), (22) is the 
following. The equation (21) guarantees that the phase 
shifts Si exactly reproduce the experimental phase shifts 
at the energies E — E\. The equation (22) fixes the 
derivatives of the phase shifts ^ at the energies E = E\ 
fitting them exactly to the derivatives of the experimental 
phase shifts at the same energies. 

The solutions E\ and (Af\X), X = 0, 1, ... , Af depend 
strongly on the values of the oscillator spacing ft£l and Af, 
the size of the inverse scattering potential matrix. Larger 
values of Af and/or ft£l, imply a larger energy interval 
< E < i?max where the phase shifts are reproduced by 
the J-matrix parameterization (7). 

A Hcrmitian Hamiltonian generates a set of normalized 
eigenvectors (n\X) fitting the completeness relation, 



£|<A0A>| 2 = 1. 



(24) 



A=0 



C$(E) - Qmn{E) Th tAr+1 C%l hl {E) - 0, (20) 



(+) 



Experimental phase shifts generate a set of {Af\\), A = 0, 
1, ... , Af that usually does not fit Eq. (24). It is likely 
that the interval of energy values used to find the sets 
of E\ and (Af\X), spreads beyond the thresholds where 
new channels are opened. Thus inelastic channels are 
present in the system suggesting the Hamiltonian should 
become non-Hcrmitian. The approach proposed in Ref. 
[4], suggests to fit Eq. (24) by changing the value of 
the component (Af\X — Af) corresponding to the largest 
among the energies E\ with X = Af . This energy E\=jy is 
usually larger than E max , the maximal energy in the in- 
terval < E < E max where the experimental phase shifts 
are available. Therefore changing (Af\X — Af) should not 
spoil the phase shift description in the desired interval of 
energies below E max ; more over, one can also vary sub- 
sequently the energy £^A=Af to improve the description of 
the phase shifts in the interval < E < E max . 

We are not discussing here the construction of the in- 
verse scattering potential but point the interested reader 
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to Ref. [4]. We note only that if the construction of 
the J-matrix inverse scattering potential is desired, one 
should definitely fit Eq. (24), otherwise the construction 
of the Hermitian interaction is impossible. In our appli- 
cations to na and pa scattering we are interested only in 
the J-matrix parameterization of scattering phase shifts; 
hence we can avoid renormalization of the component 
(A/] A = A/"). Nevertheless, we found out that this renor- 
malization improves the phase shifts description at ener- 
gies E not close to E\ values. All the results presented 
below were obtained with the help of Eq. (24). 

B. Charged particle scattering 

In the case of a charged projectile scattered by a 
charged target, the interaction between them is a super- 
position of a short-range nuclear interaction, V Nucl , and 
the Coulomb interaction, V c ': 

V = V Nucl + V c . (25) 

The Coulomb interaction between proton and nucleus is 
conventionally described as (see, e. g., [11]) 

yC ^ Ze2 erf(r/^ 
r 

In the case of pa scattering discussed below, Z — 2 and 
x = 1.64 fm [11]. 

The long-range Coulomb interaction (26) requires 
some modification of the oscillator-basis J-matrix for- 
malism described in the previous subsection. In the case 
of charged particle scattering, the wave function ui(E, r) 
at asymptotically large distances takes a form: 

ui {E, r) = k [cos Si fi ((, kr) - sin 5i g t (C, fcr)] , (27) 

where 

fi(C,kr) = ^F l (C,kr), (28) 

gi(C,kr) = -^G l ((:,kr), (29) 

Fi((, kr) and G/(C, kr) are regular and irregular Coulomb 
functions respectively, and Sommerfeld parameter ( = 
Ze 2 m/k. Instead of functions S^(E,r) and ^(E, r), one 
can introduce functions ,^(EX,r) and &(E,(,r) defin- 
ing them as infinite series, 

oo 

&{E, C,r) = J2 F m(E, Rm (r) =kf, ((, kr) (30) 

n=0 

and 

oo 

^(E,C,r) = Y,Gni(E,C)Rni(r) — -k m ({,kr), 

(31) 



in order to use i 7 ^; (£?,£) and G n i(E,Q in constructing 
continuum spectrum wave functions by means of Eq. (2). 
Such an approach was proposed by the Kiev group in 
Ref. [12]. Within this approach, the J-matrix matching 
condition at n = Af becomes much more complicated, 
resulting in difficulties in designing an inverse scattering 
approach and in shell model applications. In practical 
calculations, the approach of Ref. [12] requires the use 
of much larger values of TV", i. e. a huge extension of the 
model space when solving the algebraic problem (12), 
that makes it incompatible with the shell model appli- 
cations. Therefore it is desirable to find another way 
to extend our approach on the case of charged particle 
scattering. 

We use here the formalism of Ref. [7] to allow for the 
Coulomb interaction in the oscillator-basis J-matrix the- 
ory. The idea of the approach is very simple. Suppose 
there are a long-range V and a short-range V sh poten- 
tials that are indistinguishable at distances < r < b. In 
this case, the potential V Sh generates a wave function fit- 
ting exactly (up to an overall normalization factor) that 
of the long-range potential V at r < b. If the only dif- 
ference between V and V Sh at distances r > b is the 
Coulomb interaction, then one can equate the logarith- 
mic derivatives of their wave functions at r — b and use 
the resulting equation to express the long-range potential 
phase shifts Si in terms of the short-range potential phase 
shifts 5f h or vice versa. Note that the phase shifts 8f h 
can be obtained within the standard J-matrix approach 
discussed in the previous subsection. The recalculation 
of the phase shifts 8f h into 5i (or vice versa) appears to 
be the only essential addition in formulating such a direct 
(or inverse) Coulomb-extended J-matrix formalism. 

To implement this idea, we introduce a channel radius 
b large enough to neglect the nuclear interaction V Nucl 
at distances r > b, i. e. b > RnucI, where RnucI is the 
range of the potential V Nucl . In the asymptotic region 
r > b, the radial wave function ui(E,r) is given by Eq. 
(27). 

At short distances r < b, the wave function ui(E,r) 
coincides with uf h (E,r), the one generated by the aux- 
iliary potential 

v Sh_jV= V Nucl + V c , r<b . ( 
V ~ ) o, r > b ' - Nucl ^ ' 

obtained by truncating the Coulomb potential V c at 
r = b. The wave function uf h {E 1 r) behaves asymp- 
totically as a wave function obtained with a short-range 
interaction, 

uf h (E,r) = k[cosS^ h ji(kr)-sm8f h ni(kr)}, b > R Nucl . 

(33) 

The J-matrix formalism described in the previous subsec- 
tion, should be used to calculate the function uf h (E, r), 
the auxiliary phase shift Sf h and the respective auxiliary 
S- matrix S sh . 

Matching the functions m (E, r) and uf h (E, r) at r = b, 
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the phase shift Si can be expressed through S^ h [7]: 
+fl c _ W h { 3 i,fi)-W b {n h fi)t a nSf h 

where quasi- Wronskian 



-ji(kr) —fi((,kr) 



(35) 



r—b 



and W b {ni,fi), W b (ji,gi) and W b (ni,gi) are expressed 
similarly. The 5-matrix is given by 



Wb(v>sr) 



W 6 (fc+ gf)S Sh 



W b (hi,g+)-W b (h+,g+)S^ 



(36) 



where hf(kr) = — nj(fcr) ± iji(kr), gf{C,,kr) = 
—gi((,kr) ± ifi((,kr), and the quasi- Wronskians 
W&(/ij t ) fl'j t ) are defined by analogy with Eq. (35). The 
S-matrix poles are obtained by solving the equation 

5 +) - W b (h+,g+)S Sh = (37) 

in the complex energy plane. 

This formalism involves a free parameter, the channel 
radius b, used for construction of the auxiliary potential 
V sh . As mentioned above, b should be taken larger than 
the range of the short-range nuclear interaction y Nucl . 
On the other hand, the truncated (W + 1) x (TV + 1) 
Hamiltonian matrix H nn , (n,n' = 0, 1, ... , TV) used 
to calculate the sets of eigenvalues E\ and eigenvectors 
(A/] A) by solving the algebraic problem (12), should carry 
information about the jump of potential V sh at the point 
r = b. Therefore b should be chosen less than approx- 
imately rj\/, the classical turning point of the oscillator 
function Rj^i(r), the function with the largest range in 
the set of oscillator functions R n i (r) , n = 0, 1 , . .. , Af used 
for the construction of the truncated Hamiltonian matrix 
H l nn , (n,n' < AT). In a practical calculation, one should 
study convergence with a set of b values and pick up the 
b value providing the most stable and best-converged re- 
sults. As shown in Ref. [7], the phase shift Si calculated 
at some energy E as a function of channel radius 6, usu- 
ally has a plateau in the interval RnucI < b < rfy that 
reproduces well the exact values of Si. 

In the inverse scattering approach, first, we fix a value 
of the channel radius b and transform experimental phase 
shifts Si into the set of auxiliary phase shifts 5f h : 

, s sh _ W b (jiJi)-W b (ji,gi)tanSi 

tan<>i ~ W b (niJi)-W h (ni,gi)t a n6i- W 

Equation (38) can be easily obtained by inverting Eq. 
(34). Next, we employ the inverse scattering approach 
of the previous subsection to calculate the sets of E\ 
and (A/"|A) using auxiliary phase shifts Sf as an input. 
The J-matrix parameterization of the phase shifts Si is 
given by Eq. (34), the S- matrix poles can be calculated 
through Eq. (37). 



C. J-matrix and the shell model 

Up to this point we have been discussing the J-matrix 
formalism supposing the colliding particles to be struc- 
tureless. In applications to the na and pa scattering 
and relating the respective J-matrix inverse scattering 
results to the shell model, we should have in mind that 
the a particle consists of 4 nucleons identical to the scat- 
tered nucleon and the five-nucleon wave function should 
be antisymmetrized. The J-matrix solutions and the ex- 
pressions (7) for the phase shifts and (8) for the S'-matrix 
[or expressions (34) and (36) in the case when both the 
projectile and the target are charged] , can be used in the 
case of scattering of complex systems comprising identi- 
cal fcrmions. The components (7V"|A) entering expression 
(13) for the function Gnn{E) become, of course, much 
more complicated: they now appear to be some particu- 
lar components of the many-body eigenvector. However, 
we are not interested here in the microscopic many-body 
structure of the components (A/"|A); we shall obtain them 
by fitting the na and pa phase shifts in the J-matrix 
inverse scattering approach. 

We focus our attention here on other important ingre- 
dients entering expression (13) for Qj^jyf(E), the eigenen- 
ergies E\, related to the energies of the states in the com- 
bined many-body system, i. c. in the 5 He or 5 Li nucleus 
in the case of na or pa scattering respectively, obtained 
in the shell model or any other many-body approach uti- 
lizing the oscillator basis. One should have however in 
mind that E\ entering Eq. (13) correspond to the kinetic 
energy of relative motion, i. e. they are always positive, 
while many-body microscopic approaches generate eigen- 
states with absolute energies, e. g. all the states in 5 He 
and 5 Li with excitation energies below approximately 28 
MeV (the a-particle binding energy) will be generated 
negative. Therefore, before comparing with the set of 
E\ values, one should perform a simple recalculation of 
the shell model eigenenergies by adding to them the 4 He 
binding energy; or alternatively one can use the set of E\ 
values to calculate the respective set of energies defined 
according to the shell model definitions by subtracting 
the 4 He binding energy from each of E\. The physical 
meaning of transforming these to the shell model scale 
of values for E\ is to provide the values required from 
shell model calculations in order to reproduce the desired 
phase shifts. 

The comparison of the inverse scattering J-matrix 
analysis with the shell model results is useful, of course, 
only if the same Ml value is used both in the J-matrix 
and in the shell model and model spaces of these ap- 
proaches are properly correlated. A traditional notation 
for the model space within the shell model is N max hfl 
where N max is the excitation oscillator quanta. In the 
case of the J-matrix, we use, also traditionally, A/", the 
principal quantum number of the highest oscillator func- 
tion Rj^fi{r) included in the 'interaction region' of the 
oscillator model space where the potential energy matrix 
elements are retained. The following expressions relate 
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N max and Af in the cases of § and I partial waves (p 
waves) and ^ + partial wave (s wave): 

3 - i - 

N max = 2J\f, Af = 0,1,..., - and - partial waves, 

(39) 



N max = 2Af — 1, A/" =1,2,..., - partial wave. 



(40) 



Below we are using shell model type N max hfl notations 
for labeling both J-matrix and shell model results. 



III. ANALYSIS OF no. SCATTERING PHASE 
SHIFTS 
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A. ^ phase shifts 



We start discussion of our J-matrix analysis of na scat- 
tering with the | phase shifts. 

The J-matrix inverse scattering approach was well- 
tested in the nucleon-nucleon (NN) scattering in Ref. 
[4]. In the case of AW-scattering, the phase shifts are 
well established in a wide range of energies up to 350 
MeV in laboratory system. The goal of Ref. [4] was 
to fit scattering phase shift in the entire interval of en- 
ergies < Ei a b < 350 MeV using the smallest possible 
potential matrices or, equivalcntly, the smallest possible 
values of N max (or Af) . In the case of nucleon-a scatter- 
ing, the phase shifts are known in a small energy interval 
up to Eiab = 20 MeV and in some cases up to 25 MeV. 
On the other hand, to compare the J-matrix analysis 
with the shell model results, we are interested in large 
enough values of -/V max and in HQ values reasonable for 
shell model applications. As a result, we face a problem 
of insufficient data: some solutions E\ of Eq. (21) should 
be allowed far outside the interval of known phase shifts 
5i. The required phase shifts should be known together 
with their derivatives at the energies around E = E\ just 
to find these solutions E\ [see Eq. (14)] and respective 
eigenvector components (7V|A) [see Eqs. (22) and (23)]. 

We address the problem of data insufficiency by an 
extrapolation of the data outside the energy interval of 
known phase shifts. The J-matrix parameterizations pre- 
sented in Fig. 1 were obtained with Ml = 20 MeV in 
various model spaces. In each case two different ex- 
trapolations were used for the phase shifts at energies 
Ei a b > 20 MeV, however the experimental phase shifts 
below Ei a b = 20 MeV are equivalently well described if 
Amax is large enough. The deviation of the parameteri- 
zation from the experiment is seen at energies Ei a \> > 10 
MeV only in the case of the 2M7 model space, the small- 
est among all model spaces presented in Fig. 1, and even 
in this case the deviation is small enough. This is not 
surprising since the phase shifts given by Eqs. (7) and 
(13) in the low-energy interval are governed mostly by 
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FIG. 1: (Color online) Experimental | na phase shifts from 
Refs. [13] (stars) and [14] (filled squares) and J-matrix pa- 
rameterizations with two arbitrary extrapolations of phase 
shifts for E lab > 20 MeV obtained with hfl = 20 MeV in 
2hil, 4hQ and 6h£l model spaces. Both panels present the 
same results in different energy scales. 



the E\ values from the same interval and by the respec- 
tive eigenvector components (A/"|A). These E\ and (AA|A) 
values are determined by Eqs. (21) and (22) locally, i. e. 
they are independent from the phase shift extrapolation. 
Note that in the case of the 2MI model space, both E\ 
values lie in the energy interval of known phases, hence 
the parameterization in this model space is completely 
independent from the extrapolation and the two param- 
eterizations obtained in this model space coincide. 

The resonance energy E res and width T calculated by 
locating the S'-matrix pole by solving Eq. (20), are seen 
from Table I to be very stable and insensitive to the ex- 
trapolation of the phase shifts. 

The same insensitivity of the phase description in the 
desired energy interval to the phase shift extension out- 
side this interval, is also inherent for other na partial 
waves. Hence, we shell not waste space by discussing 
this issue in the respective subsections below. We should 
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TABLE I: The energy E res and width F (both in MeV) of 
the | resonance in the na scattering obtained with HQ = 
20 MeV in various model spaces -/V max /ifi with two different 
extrapolations of the phase shifts. 





Extrapolat 


ion 1 


Extrapolat 


ion 2 


■^V max 


E res 


r 


E res 


r 


6 


0.7713 


0.6437 


0.7718 


0.6435 


8 


0.7719 


0.6451 


0.7715 


0.6454 


10 


0.7707 


0.6417 


0.7708 


0.6416 



probably just note here that the need for data extrap- 
olation arose only due to our desire to compare the J- 
matrix results with the shell model ones; it is this desire 
that pushes us to use large enough model spaces and Ml 
values. If one is interested only in getting a high qual- 
ity J-matrix parameterization of the phase shifts and in 
extracting resonance parameters, smaller model spaces 
and/or smaller Ml values can be used without a loss of 
accuracy and without a need to have phase shifts outside 
the experimentally known energy interval. 

A resulting practical approach is to extrapolate the 
phase shifts in any reasonable way outside the energy 
interval where they are known in order to obtain the J- 
matrix parameterization of the phase shifts within this 
interval of known phases and to derive resonance param- 
eters in the same energy interval. Using such extrapola- 
tion, we study a dependence of the J-matrix phase shift 
parameterization on the size of the model space. As is 
seen from Fig. 2, larger model spaces make it possible to 
describe the extrapolated phase shifts up to larger ener- 
gies. The experimental data are perfectly reproduced in 
AMI and larger model spaces. We fail to reproduce the 
experiment for Ei a b > 12 MeV in the 2MI model space. 
Note however that deviations from the experiment are 
not very large and we obtain a very good description of 
the phase shifts at laboratory energies below 12 MeV in- 
cluding the resonance region. The smallest possible 0M7 
model space fails to provide a reasonable description of 
the phase shifts at all energies. 

As mentioned above, the description of the phase shifts 
can be extended to larger energies not only by using 
larger model spaces but also by using larger Ml values. 
This is illustrated by Fig. 3. Even with Ml = 5 MeV 
we manage to describe the phase shifts in the 6MI model 
space up to approximately Ei a b = 17 MeV. The descrip- 
tion of all experimentally known phase shifts is perfect 
in this model space with Ml = 10 MeV and larger. 

The results of calculations of the ^-matrix pole posi- 
tion are presented in Fig. 4. The calculated resonance 
energy E res and width T are seen to be very stable in a 
wide range of Ml values and model spaces (note a very 
detailed energy scale in Fig. 4). Our results are in a 
very good correspondence with the results of a detailed 
study of Ref. [15]. The authors of this paper performed 
Resonating Group Method calculations of Na scattering 
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FIG. 2: (Color online) The J-matrix parameterization of the 
| na phase shifts obtained with Ml = 20 MeV in various 
model spaces. Different panels present the same results in 
different energy scales. Experimental phase shifts: stars — 
Ref. [13], filled squares — Ref. [14]. 

with phenomenological Minnesota NN interaction fitted 
to reproduce with high precision the na and pa phase 
shifts and calculated the position of the S'-matrix pole. 
The extended multichannel i?-matrix analysis of 5 He and 
5 Li including two-body channels N+a and d+t or d+ 3 He 
along with pseudo-two-body configurations to represent 
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FIG. 3: (Color online) The J-matrix parameterization of the § na phase shifts obtained in the QhQ model space with different 
values of oscillator spacing hQ. See Fig. 2 for details. 



0.80 - - 



0.79- 



> 

0.78 
0.77 



0.76 -- 



Extended R-matrix 

— T — 6hn 
-A- I2hn 



S-matrix, RGM 

—\ I I I I I I I I I I I I r~ 



0.80 -- 



0.79 



> 



=- 0.78 
ilT 



0.77 



0.76 



10 15 

hQ (MeV) 
(c) 



20 



Extended R-matrix 



-A- hn = 10MeV 
— ▼ — hn = 20 MeV 




0.66 



0.65 



0.64 



0.63 - - 



(b) 



Extended R-matrix 



X 



y 



— T— 6hn 
-A- I2hn 



S-matrix, RGM 




8 10 12 



FIG. 4: (Color online) The na | resonance energy in the center-of-mass frame (left) and width (right) obtained by calculating 
the position of the S-matrix pole by means of the J-matrix parameterizations with different h£l values (upper panels) and in 
different model spaces (lower panels) . Horizontal lines present the results of Ref. [15] : the analysis of the resonance parameters 
in the extended _R-matrix approach (dashed) and calculations of the S-matrix pole position in the Resonating Group Method 
(dash-dotted). 
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FIG. 5: (Color online) The lowest state i?A=o obtained in the 
J-matrix parameterization of the § na phase shifts with 
different Ml values in various model spaces. 

the breakup channels n+p+t or n+p+ 3 He, was also per- 
formed in Ref. [15] using data of various authors on the 
differential elastic scattering cross sections, polarization, 
analyzing-power and polarization-transfer measurements 
together with neutron total cross sections. Our very sim- 
ple J-matrix analysis utilizing only the elastic scattering 
phase shifts, is competitive in quality of resonance pa- 
rameter description with these extended studies of Ref. 
[15]. 

We note that while the phase shifts and resonance pa- 
rameters are very stable, the energies E\ entering Eq. 
(13) vary essentially with Ml and model space. In par- 
ticular, this is true for the lowest of these energies E\ =0 
shown in Fig. 5 (note a very large difference in energy 
scales in Figs. 4 and 5). This energy being obtained 
in shell model studies, would be associated traditionally 
with the resonance energy E res . Such a conventional 
association is clearly incorrect: this lowest eigenstate 
i?A=o differs significantly in energy from E res while the 
phase shifts and resonance energy and width are well 
reproduced; just this energy E\ = q, very different from 
E res , is needed to have a perfect description of scatter- 
ing data and resonance parameters including E res itself. 
The E\=o dependencies of the type shown in Fig. 5 are 
inherent in other partial waves and in the case of pa 
scattering. We study the E\ =Q dependencies on Ml and 
model space in more detail below in Section V where we 
compare them with the results of No-core Shell Model 
calculations. 




E , ab ( MeV ) 
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FIG. 6: (Color online) The J-matrix parameterization of the 
i na phase shifts obtained with Ml — 20 MeV in various 
model spaces. See Fig. 2 for details. 



B. ^ phase shifts 

We present in Figs. 6 and 7 respectively the J-matrix 
parameterizations of na h phase shifts obtained with 
the same Ml in different model spaces and with different 
Ml values in the same model space. The description of 



the 2 phase shifts with different Ml values and in dif- 
ferent model spaces follows the same patterns as in the 
case of the | phase shifts. The only difference is that 
a high-quality description of the phase shifts at energies 
Eiab > 10 MeV is attained in larger model spaces. How- 
ever, in the 8htt and larger model spaces the description 
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FIG. 7: (Color online) The J-matrix parameterization of the ^ na phase shifts obtained in the (SKI model space with different 
values of oscillator spacing Ml. See Fig. 3 for details. 
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FIG. 8: (Color online) The na | resonance energy in the center-of-mass frame (left) and width (right) obtained by calculating 
the position of the S-matrix pole by means of the J-matrix parameterizations with different Ml values (upper panels) and in 
different model spaces (lower panels). See Fig. 4 for details. 
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of all known phase shifts is perfect. 

Figure 8 presents the results of our calculations of the 
i resonance energy and width. The variations of E res 
and r with increasing Ml or model space are larger than 
in the case of the | resonance; note however that the 

energy of the | resonance and its width are also much 
larger. At any rate, the variations of resonance param- 
eters are not large and our results for E res and T are 
stable enough with respect to the choice of Ml value and 
model space. The energy and width of the ^ resonance 
also compare well with the results of Ref. [15]. 



C. i phase shifts 

In describing the i + phase shifts, one should have in 
mind that the lowest s states are occupied in the a- 
particle and due to the Pauli principle these states should 
be inaccessible to the scattered nucleon. There are two 
conventional approaches to the problem of the Pauli for- 
bidden s state in the n + a system. The first approach is 
to add a phenomenological repulsive term to the s wave 
component of the na potential (see, e. g., [16]). This 
phenomenological repulsion excludes the Pauli forbidden 
state in the n + a system and is supposed to simulate 
the Pauli principle effects in more complicated cluster 
systems. Another approach is to use deep attractive na 
potentials that support the Pauli forbidden s state in the 
n + a system (see [11, 17, 18]). In the cluster model stud- 
ies, the Pauli forbidden state is excluded by projecting it 
out [11, 18, 19]. 

In our J-matrix inverse scattering approach, we can 
simulate both the potentials with repulsive core and with 
a forbidden state. In the first case, when the system does 
not have a bound state, we go on with the same proce- 
dure as in the above cases of 1 and | partial waves; 

the energy dependence of the input \ + phase shifts is 
responsible for generating proper details of the na inter- 
action potential matrix. In the other case, the simplest 
way to simulate the presence of the forbidden state in 
the system is to suppose that this state is described by 
a pure Osj/a oscillator wave function. The energy of the 
forbidden state is equal in this case to the Hamiltonian 
matrix element Hqq° which is of no interest for us in this 
study, all the matrix elements Hj^° and should be 

set equal to zero to guarantee the orthogonality of the 
forbidden state to scattering states which have the wave 
functions given by the expansion (2) where the Osi/2 os- 
cillator state is missing, i. e. a n= oj=o{E) = for all 
energies E > 0. Within this model, the forbidden state 
[20] does not contribute to the function GaTm(E) [see Eq. 
(13)] since the component (Af\ \ = 0) = 0. In the inverse 
scattering approach, we use the first J\f solutions of Eq. 
(21) disregarding the highest in energy solution 
while constructing the function Gnn(E). 

In Fig. 9 we present the J-matrix parameterization of 
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FIG. 9: (Color online) The J-matrix parameterization of the 
na phase shifts obtained in the 7 Ml model space with 
different values of oscillator spacing Ml. See Fig. 3 for details. 

the | + phase shifts in elastic na scattering in the 7 Ml 
model space with different values of the oscillator spac- 
ing Ml. As usual, larger Ml value makes it possible to 
describe the phase shifts in a larger energy interval. A 
new and interesting issue is the difference in behavior of 
the phase shifts in the models with and without a forbid- 
den state. A more realistic model with forbidden state 
provides a proper dependence of the phase shifts: start- 
ing with 180° at zero energy, they tend to zero at large 
energies. The forbidden state makes the same contribu- 
tion to the Levinson theorem as any other bound state 
providing the 180° difference between the phase shifts at 
zero and infinite energies. The model without a forbid- 
den state generates the phase shifts returning at large 
energies back to their zero energy value. In what follows, 
we use the potential model with a forbidden state. Note 
however that in the energy interval of known phase shifts, 
the parameterizations of both models are indistinguish- 
able. The E\ values provided by both models in this 
energy interval, are the same. 
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FIG. 10: (Color online) The J-matrix parameterization of the 
j + net phase shifts obtained in the model with forbidden state 
with h£l — 20 MeV in various model spaces. See Fig. 3 for 
details. 

The j + phase shifts parameterizations in different 
model spaces with Ml = 20 MeV perfectly describe the 
data (Fig. 10). At larger energies, they follow general 
trends: smaller model spaces result in a faster fall off of 
the phase shifts to zero value. 



IV. ANALYSIS OF pa SCATTERING PHASE 
SHIFTS 

A. | phase shifts 

The J-matrix approach to pa scattering involves an 
additional parameter 6, the channel radius used to define 
the auxiliary potential V Sh by truncating the Coulomb 
interaction at r = b [see Eq. (32)]. We start our discus- 
sion of the J-matrix inverse scattering description of pa 
scattering from the analysis of the 6-dependence of the 
| pa phase shift parameterization. 
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FIG. 11: (Color online) The J-matrix parameterization of the 
| pa phase shifts obtained with ftfi = 20 MeV in the lOftfi 
model space with various b values. Two panels present the 
same results in different scales. Experimental phase shifts: 
open squares — Ref. [21], open circles — Ref. [22], crosses — 
Ref. [23]. 



We present in Fig. 11 the pa | phase shift parame- 
terizations obtained with different channel radii b in the 
lOhfl model space with HVL = 20 MeV. The experimen- 
tal data are seen to be perfectly described in the interval 
of b values 6 fm < b < 10 fm. However we did not 
find a way to reproduce accurately the phase shifts with 
b > 11 fm, in particular at energies between 10 and 20 
MeV. This is not surprising since the classical turning 
point of the highest oscillator function Rj^i{r) involved 
in the construction of the truncated Hamiltonian H nn , 

(n, n' < AT), rff = 8.06 fm in this case. The | reso- 
nance energy and width dependences on b obtained from 
these parameterizations, are shown in Fig. 12. The res- 
onance parameters are seen to be stable enough with b 
varying between 6 and 10 fm. The b dependence of the 
energy E\=o of the lowest state obtained in the J-matrix 
parameterization, has a plateau between 7 and 10 fm (see 
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FIG. 12: (Color online) Dependence of the | pa resonance 
energy in the center-of-mass frame (upper panel) and width 
(lower panel) on the channel radius b in the lOftO model space 
calculations with hu> = 20 MeV. See Fig. 4 for details. 



Fig. 13). 

The bottom line of these studies is that the results 
are nearly 6-independent for b values in some vicinity of 
the classical turning point r&. This conclusion remains 
valid for other partial waves of the pa scattering and we 
are not discussing 6-dependences in the following subsec- 
tions. The remainder of the calculations presented here 
are performed with b = rfy. 

We present in Fig. 14 the J-matrix parameterization of 
the | pa phase shifts obtained in various model spaces 
with hil — 20 MeV. The data are well-described in 4M7 
and higher model spaces. Some deviation from experi- 
ment is seen only for the 2MI model space starting from 
laboratory energies about 20 MeV. However, the reso- 
nance region is perfectly described even in this very small 
2/iSl model space as is seen from the lower panel of Fig. 
14 where the enlarged energy scale is used. The J-matrix 
parameterization is also insensitive to the variation of the 
HQ value in the whole interval of known phase shifts in- 
cluding the resonance region (see Fig. 15). Therefore it 
is not surprising that we obtain a very stable description 
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FIG. 13: The b dependence of the lowest state £a=o m the 
10hQ model space J-matrix parameterizations with hw = 20 
MeV. 
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FIG. 14: (Color online) The J-matrix parameterization of the 
pa | phase shifts obtained with h£l = 20 MeV in various 
model spaces. See Fig. 11 for details. 



14 



m 

a 

■D 



120 n 



100- 



(a) 



80- 



60- 



40- 



20- 




hQ = 10 MeV 
ha = 15 MeV 
hn = 20 MeV 
ha = 25 MeV 
ha = 30 MeV 



10 15 
E lab (MeV) 



20 



25 



120 



100 



«r 80 -i 

£ 

g 1 60-| 



40- 



20- 



(b) 




ha = 10 MeV 
ha = 15 MeV 
ha = 20 MeV 
ha = 25 MeV 
ha = 30 MeV 



E lab (MeV) 



FIG. 15: (Color online) The J-matrix parameterization of the pa | phase shifts obtained in the WhH, model space with 
various fitt values. See Fig. 11 for details. 
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FIG. 16: (Color online) The pa | resonance energy in the center-of-mass frame (left) and width (right) obtained by calculating 
the position of the S-matrix pole by means of the J-matrix parameterizations with different Ml values (upper panels) and in 
different model spaces (lower panels). See Fig. 4 for details. 
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FIG. 17: (Color online) The J-matrix parameterization of the 
pa i phase shifts obtained with Ml = 20 MeV in various 
model spaces. See Fig. 11 for details. 
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FIG. 18: (Color online) The J-matrix parameterization of the 
pa i phase shifts obtained in the 10M7 model space with 
various hQ values. See Fig. 11 for details. 

of the resonance energy and width (see Fig. 16), one that 
is independent of the model space and Ml value. 

Our results for the | resonance parameters are very 
close to the ones obtained in the analysis of Ref. [15]. 



resonance region is perfectly described. Our results for 
the resonance energy and width correspond well to the 
analysis of Ref. [15]. The resonance parameters are sta- 
ble with respect to variations of the model space and Ml 
(see Fig. 19). Of course, the variations of E res and T 
in Fig. 19 are much larger than in the case of the | 
resonance, but the \ resonance energy and width are 

also much larger than the energy and width of the | 
resonance. 



C. ^ + phase shifts 

In the case of s wave of pa scattering, we can also 
use interaction models with and without a forbidden 
state. The main features of the J-matrix parameteriza- 
tions within these models in the case of the pa scattering 
are the same as in the case of na scattering; in particu- 
lar, the phase shift description in the low-energy region 
covering the whole region of known phase shifts, is iden- 
tical within these interaction models. In what follows, 
we present only the results obtained in the model with 
forbidden state which we suppose to be more realistic. 

The J-matrix parameterizations of the pa ^ + phase 
shifts obtained in various model spaces with Ml = 20 
MeV, are presented in Fig. 20 in two scales. The low- 
energy phase shifts up to approximately Ei a b — 10 MeV 
are perfectly reproduced in all model spaces. Starting 
from Ei a b = 10 MeV, there are some deviations from the 
experiment that are well seen in the right panel of Fig. 20 
where a larger scale is used. Surprisingly, the deviations 
from experimental phase shifts are larger in larger model 
spaces. The deviations are not large but not negligible. 

The J-matrix parameterizations obtained with various 
Ml values in the 11MI model space, are shown in Fig. 
21. The theoretical curves are nearly indistinguishable 
below Eiab = 10 MeV reproducing well the experimental 
data. Some difference between parameterizations is seen 
in the high-energy part of the interval of known phase 
shifts. All J-matrix parameterizations presented in Fig. 
21 reasonably describe the phenomenological data in the 
whole energy interval of known phase shifts. The worst 
description of the phase shifts in the 11M7 model space 
is obtained with Ml = 15 MeV. 



B. ^ phase shifts 

We obtain a high-quality J-matrix parameterization of 
the pa ^ phase shifts, very stable with variations of the 
model space or oscillator spacing Ml. A small deviation 
from the experiment at large energies is seen in Fig. 17 
in the 2MI model space only. The parameterizations ob- 
tained in the 1QMI model space with Ml values ranging 
from 10 to 30 MeV, are indistinguishable in Fig. 18. The 



V. J-MATRIX AND SHELL MODEL 
EIGENSTATES 

Up to now, we were discussing the J-matrix inverse 
scattering description of scattering observables in the 
n + a and p + a nuclear systems. It is very interesting 
to investigate whether these observables correlate with 
the shell model predictions for 5 He and 5 Li nuclei. It 
should be done, as we have shown above, by comparing 
the eigenenergies E\ obtained in the J-matrix inverse 
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FIG. 19: (Color online) The pa | resonance energy in the center-of-mass frame (left) and width (right) obtained by calculating 
the position of the S-matrix pole by means of the J-matrix parameterizations with different Ml values (upper panels) and in 
different model spaces (lower panels). See Fig. 4 for details. 
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FIG. 20: (Color online) The J-matrix parameterization of the pa | + phase shifts obtained in the model with forbidden state 
with hQ. = 20 MeV in various model spaces. See Fig. 11 for details. 
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FIG. 21: (Color online) The J-matrix parameterization of the 
pa | + phase shifts obtained in the model with forbidden state 
in the llftfi model space with various hQ values. See Fig. 11 
for details. 



scattering approach with the energies of the states ob- 
tained in the shell model. 

We calculate the lowest 5 He and 5 Li states of a given 
spin and parity in the No-core Shell Model approach 
[8] using the code MFDn [24] and the JISP16 nucleon- 
nucleon interaction [6, 25]. We do not make use of ef- 
fective interactions calculated within Lee-Suzuki or any 
other approach. That is, all results presented here are 
obtained with the 'bare' JISP16 NN interaction which 
is known [6, 29, 30] to provide a reasonable convergence 
as basis space size increases. One may note that the 
No-Core Shell Model with a bare interaction and with a 
truncated configuration basis may also be referred to as a 
"configuration interaction" or "CI" type calculation [31]. 

In all cases, the calculations of the 4 He ground state 
energy is performed with the same Ml value and in the 
same N meix hfl model space. These 4 He ground state en- 
ergies are used to calculate the reaction threshold while 
comparing the J-matrix E\ values (defined with regard 
to the reaction threshold) with the shell model results. 
Therefore our reaction threshold is model space and Ml- 
dependent, however these dependencies are strongly sup- 



pressed in large enough model spaces. This definition of 
the reaction threshold is, of course, somewhat arbitrary. 
We use it supposing that our definition provides a consis- 
tent way to generate energies relative to the 4 He ground 
state energy within the No-core Shell Model approach 
employing a finite basis. 

The No-core Shell Model results for the lowest 5 He 
and 5 Li | , i and h states are compared with the 
respective J-matrix E\ = q values in Figs. 22 and 23. For 
each spin and parity, the J-matrix £\=o values obtained 
with the same Ml value, are seen to decrease with in- 
creasing model space (see also Fig. 5); the same model 
space dependence is well-known to be inherent for the 
shell model eigenstates. However the Ml dependences of 
the J-matrix E\ = q and shell model eigenstates, are very 
different: the shell model eigenstates are known to have 
a minimum at some Ml value while the inverse scattering 
E\ =a are seen from Figs. 22 and 23 to increase nearly lin- 
early with Ml; the slope of the Ml dependence of E\ = o is 
larger for wider resonances. As a result, the shell model 
predictions differ from the results of the inverse scattering 
analysis for small enough Ml values. However, a remark- 
able correspondence between the shell model and inverse 
scattering results is seen at large enough Ml values start- 
ing from approximately Ml — 20 MeV. The agreement 
between the shell model and J-matrix inverse scattering 
analysis is improved with increasing model space; it is 
probable that this is partly due to the improvement in 
larger model spaces of the calculated threshold energy in 
our approach. The shell model description of the lowest 



and i states is somewhat better than the lowest | 



state description in both 5 He and 5 Li nuclear systems. 
The lowest | state description is however not so bad 

(note a more detailed energy scale for the § state in 
Figs. 22 and 23): the difference between the shell model 
predictions and the J-matrix analysis results is about 0.5 
MeV in large enough model spaces and for large enough 
Ml values. An excellent description of the | states in 
5 He and 5 Li combined with some deficiency in description 
of the | states in the same nuclei, is probably a signal 
of a somewhat underestimated strength of the spin-orbit 
interaction generated by the JISP16 NN interaction in 
the p shell. 

We suppose that the results presented here illustrate 
well the power of the proposed J-matrix analysis, a new 
method that makes it possible to verify a consistency of 
shell model results with experimental phase shifts. To 
the best of our knowledge, this is the only method which 
can relate the shell model results to the scattering data 
in the case of non-resonant scattering like the i + na 
and pa scattering. In the case of negative parity reso- 
nances in 5 He and 5 Li discussed here, the J-matrix anal- 
ysis generally suggests that the shell model should gen- 
erate the respective states above the resonance energies 
supplemented by their widths. Note that the J-matrix 
E\ = q only in some cases lie inside shaded areas showing 
the resonance energies together with their widths in Figs. 
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FIG. 22: (Color online) E\ = o values for not scattering ob- 
tained in the J-matrix inverse scattering approach in the 
center-of-mass frame (filled symbols) and respective lowest 
eigenstates of the He nucleus obtained in the No-core Shell 
Model (respective empty symbols). The resonance energies 
together with their widths are shown by shaded areas. 



FIG. 23: (Color online) E\ = o values for pa scattering ob- 
tained in the J-matrix inverse scattering approach in the 
center-of-mass frame and respective lowest eigenstates of the 
Li nucleus obtained in the No-core Shell Model. See Fig. 22 
for details. 
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22 and 23, and in all these cases, the intersection of the 
E\ = q with the resonance is seen only at small enough Ml 
values where the shell model predictions fail to follow the 
J-matrix analysis results. This is a clear indication that 
one should be very accurate in relating the shell model 
results to the resonance energies, at least in the case of 
wide enough resonances. 

VI. CONCLUSIONS 

We suggest a method of J-matrix inverse scattering 
analysis of elastic scattering phase shifts and test this 
method in applications to na and pa elastic scattering. 
We demonstrate that the method is able to reproduce 
| , i and i + na and pa elastic scattering phase shifts 
with high accuracy in a wide range of the parameters of 
the method like the oscillator spacing Ml, model space 
and the channel radius b in the case of pa scattering. The 
method is very simple in applications, it involves only a 
numerical solution of a simple transcendental equation 
(21). 

When the J-matrix phase shift parameterization is ob- 
tained, the resonance parameters, resonance energy and 
width, can be obtained by locating the S'-matrix pole by 
solving numerically another simple transcendental equa- 
tion (20). The resonance energies and widths are shown 
to be stable when Ml or other J-matrix parameters are 
varied. Our results for | and ^ resonant states in 5 He 
and 5 Li are compared in Table II with the results of other 
authors. Our results are in line with the results of other 
studies; in general, the better agreement is seen with Rcf. 
[15], the most recent among all publications presented in 
the Table. Csoto and Hale performed two different anal- 
yses in Ref. [15]: (i) RGM search for the S'-matrix poles 
based on a complicated enough calculations within the 
Resonating Group Model with effective Minnesota NN 
interaction fitted to the nucleon-a phase shifts, and (ii) 
Extended i?-matrix analysis of 5 He and 5 Li including not 
only N + a channel but also d + t or d + 3 He channels 
along with pseudo-two-body configurations to represent 
the breakup channels n + p + t or n + p + 3 He and us- 
ing a wide range of data on various reactions. We note 
that our very simple J-matrix approach uses only a very 
limited set of data as an input, na or pa phase shifts. 
We suppose that the proposed approach can be useful 
in analysis of elastic scattering in other nuclear systems 
and serve as an alternative to the conventional i?-matrix 
analysis. 

A very interesting and important output of the J- 
matrix inverse scattering analysis of the phase shifts is 
the set of E\ values which are directly related to the 
eigenenergies obtained in the shell model or any other 
model utilizing the oscillator basis, for example, the Res- 



onating Group Model. The J-matrix paramctcrizations 
provide the energies of the states that should be obtained 
in the shell model or Resonating Group Model to gener- 
ate the given phase shifts. These energies are shown to be 
model space and M7-dependent and very different from 
the energies of at least wide enough resonances which are 
conventionally used to compare with the shell model re- 
sults. More, the J-matrix analysis is shown to provide 
the shell model energies even in the case of non-resonant 
scattering such as the i + nucleon-a scattering. 

Our comparison of the lowest E\ = o with the No-core 
Shell Model results shows that the shell model fails to re- 
produce the phase shifts if small Ml values are employed 
in the calculations. When Ml and/or model space size 
is increased, the shell model predictions approach E\ = q 
values obtained in the J-matrix signaling that the shell 
model results become more and more consistent with the 
experimental phase shifts. However some difference be- 
tween the No-core Shell Model predictions and the J- 
matrix analysis results is seen even in the largest model 
spaces used in this study. This difference is really not 
large, its possible sources are the following, (i) There is 
an ambiguity in the threshold energies used to relate the 
absolute negative energies obtained in the shell model 
and positive E\ values defined relative to the reaction 
threshold, (ii) Unfortunately, there is no NN interaction 
providing correct energies for, at least, light nuclei. The 
JISP16 NN interaction is good enough and provides re- 
liable predictions for energies of levels in all s and p shell 
nuclei [6, 29, 30]. However, there are small differences 
between JISP16 level energy predictions and experiment; 
these differences are of the same order as the differences 
between the J-matrix E\ =0 values and our No-core Shell 
Model results. Probably we shall use the J-matrix re- 
sults discussed above while attempting to design a new 
improved version of the JISP16 interaction by trying to 
eliminate the discrepancy between the shell model results 
and the J-matrix analysis of nucleon-a scattering. 

Of course, the J-matrix can be used to relate the 
shell model energies and data on nucleon scattering 
by other nuclei. Generally, one can also use other 
elastic scattering data, for example, nucleus- nucleus 
clastic scattering phase shifts to get the E\ values that 
should be obtained in the shell model studies of the 
respective compound nuclear systems: the shell model 
must generate the states with the same energies in the 
same model space and with the same Ml value to have a 
chance to generate the experimental phase shifts. 
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